WATER FRACTIONS IN EXTRASOLAR PLANETESIMALS 
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With the goal of using externany-polluted white dwarfs to investigate the water 
fractions of extrasolar planetesimals, we assemble from the literature a sample that we 
estimate to be more than 60% complete of DB white dwarfs warmer than 13,000 K, 
more luminous than 3 x 10~^ Lq and within 80 pc of the Sun. When considering all 
the stars together, we find the summed mass accretion rate of heavy atoms exceeds 
that of hydrogen by over a factor of 1000. If so, this sub-population of extrasolar 
asteroids treated as an ensemble has little water and is at least a factor of 20 drier than 
CI chondrites, the most primitive meteorites. In contrast, while an apparent "excess" 
of oxygen in a single DB can be interpreted as evidence that the accreted material 
originated in a water-rich parent body, we show that at least in some cases, there can 
be sufficient uncertainties in the time history of the accretion rate that such an argument 
may be ambiguous. Regardless of the difficulty associated with interpreting the results 
from an individual object, our analysis of the population of polluted DBs provides 
indirect observational support for the theoretical view that a snow line is important in 
disks where rocky planetesimals form. 

Subject headings: planetary systems - stars, white dwarf 
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1. INTRODUCTION 



Evidence is now strong that most e xternally-polluted white dwarfs derive their hea.vy atom by 
accretion from disrupted plane tesimals ( Jura |[2003 : Gaensicke et al. Ibood : I Jura Iboosl : IParihi et al 



201C: Zuckerman et al. 



ual objects (IKleiri et al 



Mehs et al. 



2011 



20111: 



201C) 



2010 



Observational progress is rapid; both deta ile d studies of Individ 
Dufour et al. IboidJFarihi et al. Iboila 



2011 



Zuckerman et al. 



Zuckerman et al. 



2011 



and broader surveys (jKoester et al 



Vennes et al. II201C: 



2011 



Kawka et al. 



l2oiol ;i are now being performed. Here, we combine recent st udies of white 



dwar fs to assess the water fraction - either in pure form or in hydrated minerals (jRivkin et al. 
2002) - within an ensemble of extrasolar minor planets. 



^Department of Physics and Astronomy, University of California, Los Angeles CA 90095-1562; jura, 
sxu@astro.ucla.edu 

^We use "heavy atoms" to mean all elements with atomic number greater than two. We do not follow astronomical 
convention and use "metals" for this matter because of this word's different meaning in the context of the physics of 
planets. 
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Oxygen is so cosmically abundant that H2O is predicted to be a major constituent in solid 
planets that form in regions where the temperature in the planet- forming disk is sufficiently low 



that ice is stable. Water is widespread in the outer solar system (jjewitt et al. 1 120071 : lEncrenaz 



20081 ) . and can be more than 50% of the mass of Kuiper Belt Objects and comets. C eres, the 
largest asteroid, is probably ^^25 % water (IMcCord k, Sotin I l2005l : iThomas et al. I l2005l ) and CI 
chondrites are nearly 20% water (IWasson &: Kallemeyn Ml988l ). In contrast, much of the inner 
solar system is dry. For ex ample, the refractory-rich CV chondrites are only about 2% water 

Although the amount of internal water is not wel l known, the 

The usual 



(jWasson fc Kallemevn I llQSSD . 



20071 ) 



water fraction of bulk Earth is between 0.06% and 2% (jvan Thienen et al 
theoretical explanation for the large radial variation in water content in the solar system is that 
there is a snow line in the planet- forming disk and that ice forms and condenses into planetesimals 



only in the cold, outer regio i is (lEncrenaz 



planets (jPortnev et al. 



2007 



200; 



Water may be c ommon in some extr asolar 
Seaeer et al. I l2007l : ISotin et al. I I2OO7I : IValencia et all 12007 ')■ but 



actual evideiice of ice-rich extrasolar p lanets is sparse and uncertain (jCharbonneau et al. 



Gould et al. 


2010: 


Tinnev et al. 


2011 


Because most 


internal water 


can 



200s 



km d uring a 3 Mq star's pre-white-dwarf Asymptotic Giant Branch (AGB) evolution (jjura &: Xu 



2OIOI ) , measurements of oxygen in an externally-polluted whit e dwarf's atrnqsphe rqil can serve as a 



tool to assess the amount of water in an accreted parent body. iKlein et al. I (|2010l ) showed that the 
oxygen accreted by GD 40 could have been carried in oxides such as MgO and Si02 and that water 
is less than 10% by mass of th e photospheric pol lution. Similarly, in HS 2253+8023 less than 10% 
of the accreted mass is water (jKlein et al. 1120111 ). However, because of the inherent uncertainties 
in abundance determinations, there is little prospect of using this technique to reduce much below 
~10% by mass the bound on the amount of water carried in the parent bodies. 



Water-ri c h pare nt bodies might be identified if there is an excess of oxygen (jJura &: Xu II2OIOI ). 
Farihi et al. I (|2011al ) proposed that GD 61 has accreted a water-rich asteroid, but, as discussed in 
Section 3 below, we argue that the evidence in support of this suggestion is inconclusive because 
their argument is somewhat sensitive to the unknown time history of the accretion rate. 

Another pathway to study the water content in extrasolar asteroids is to measure the amount 
of hydrogen that results from p lanetesimal acc retion onto DB white dwarf^ because this light 
gas never gravitationally settles (IKoester II2009I ). At least 1/3 of DBs warmer than 14,000 K ar e 
externally-polluted by heavy atoms derived from ancient planetesimals (jZuckerman et al. II2OIOI ). 
These same stars also acquire whatever hydrogen that was bound into water in the parent body. 



Gaensicke et al. 1 2010| ) reported dredged-up oxygen in two massive, cool white dwarfs, but there is no known 



way that the systems we discuss here have internally-polluted atmospheres. 

^DBs are white dwarfs where the helium lines dominate the spectrum. Related objects are the DBAs where 
hydrogen also is detected and the DZs which are thought to have helium-dominated atmospheres but are too cool for 
helium lines to appear in the spectrum. 
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By assessing the relative amounts of hydrogen and heavy-atoms in the stellar photospheres, we 



can co nstrain the amount of water in the accreted planetesimals. For example, IZuckerman et al 



([20101) showed that the upper limit to the hydrogen abundance in G241-6 indicates that little water 
has been accreted. 

Because hydrogen accumulates over the entir e cooling age of the st ar and because the white 
dwarf may have substantial primordial hydrogen ([Bergeron et al. Il201ll ). we cannot attribute the 
observed hydrogen to any particular accretion event. Instead, we consider a well-defined ensemble 
of white dwarfs and sum the individual hydrogen accumulation rates to compare with the summed 
heavy atom accretion rates. While particular white dwar fs such as GD 362 have acquired so much 
hydrogen that they may have accreted ice-rich matter ([Jura et al. 1 12009| ) , here we consider the 
warm DBs in aggregate. By considering the ensemble as a whole, this procedure can probe much 
smaller water fractions than allowed by assessing the amount of hydrogen in individual objects. 
Below, we show that water is less than 1% of the aggregate mass in our sample of extrasolar 
asteroids, an order of magnitude improvement over what can be achieved by examining each star 
on a case-by-case basis. 

In Section 2, we describe a volume-limited sample of DBs and compare the summed accretion 
rates of both heavy atoms and hydrogen. In aggregate, the accreted material is dry. In Section 3, 
we describe models for a quasi-steady state accretion and argue that the evidence that GD 61 is 
accreting ice-rich material is inconclusive. In Section 4 we discuss some further implications of our 
analysis, and in Section 5 we present our conclusions. 



2. AGGREGATE ACCRETION RATES 

2.1. The Sample of Stars 

In order to assess the relative amounts of accreted hydrogen and heavy elements, we assemble 
a volume-limited sample of DBs that is as bias-free as possible. We cannot use the available survey 
of white dwarfs withi n 20 pc of the Sun that is nearly complete because it contains only one DB 
([Holberg et al. II2008I ). As a compromise between having enough stars to perform an analysis yet 
not being too incomplete, we adopt the outer boundary of our sample at 80 pc. Because it is 
difficult to distinguish helium in atmospheres cooler than 13,000 K and such cool DBs are difficult 
to characterize; we only consider stars warmer than this temperature. 

Using available literature, we list in Table 1 57 DBs with luminosities greater than 3 x 10~^ 
Lq or MboI = 11 mag, alon g with their dist a nces, temperatures, galactic latitudes, masses, and 
accretion rates. According to [Bergeron et al. I ( 201ll ). the space density of such DBs is 5.15 x 10~^ 



stars pc ^, and we might expect to identify 110 appropriate objects within 80 pc. We now assess 
the incompleteness of our sample in more detail. 



If the stars are distributed isotropically in the sky, we would expect 50% to lie at \b\ < 30°. 
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However, we only identify 17 stars in this zone. The incompleteness of Tab le 1 is at least partly 



(2011 



because some stars near the Galactic Plane are missed. As demonstrated bv lLepine et al 
and as ill ustrated by the recen t discovery of a DAZ at a distance of 55 pc with b = -8° and Teff = 
20,900 K (jVennes et al. Il2010l ). there are nearby warm white dwarfs yet to be found in this region 



In actuality, we expect fewer than half the stars lie at |6| > 30° because the space density of 
stars decreases away from the Galactic Plane with a scale height, h, which depends upon the mass 
of the main-seqeunce progenitor. Assume that 



n 



no e 



-\r smb\/h 



(1) 



where no is the local density in the mid-Plane. If so, then the number of stars in a spherical volume 
of radius D centered on the Sun assumed to lie in the midplane of the Galactic Plane, N{D), can 
be found from a Taylor series expansion as: 



NiD) 



3D 1 1)2 

^ ~ 8 X ^ loT^ 



(2) 



With L) = 80 pc and h no larger than 150 pc (jGilmore Zeilik Il2000l ). we expect to find no more 
than 91 DBs. Because we identify 57 stars in Table 1, our sample is at least 60% complete. 



(2011 


) and 


Voss et al. 


(2007) 



a star is considered in bot h samples, to be as consistent as possible, we adopt the stellar parameters 
in 



Bergeron et al. I (|201ll ). For a few stars, only very limited data are available. 



One of the most important DBs for our purposes is HS 2253+8023 because it has a large 
amount of externa l pollution. The star's gravity is not well constrained in the detailed study of 
Klein et al. I (|201ll ): we adopt their most probable value. The distance to this star is derived from 
its radius, effective temperature and J-band 2MASS magnitude. We ado pt a similar approach t o 



determine distances to the other stars whose distances are not provided bv lBergeron et al. I (j2011 



The average rnass of the DBs in Table 1 is 0.67 Mq. Using the initial mass final mass relation 
of I Williams et al. I ([20091), this implies a typical main-sequence progenitor mass of 2.5 Mq. 



2.2. Upper Bounds to the Hydrogen Accretion Rates 



We now consider a DB's hydrogen budget. For the z'th star, we compute the upper bound to 
the hydrogen accretion rate averaged over the entire cooling age of the star, Mi(H), simply as the 
total mass of accumulated hydrogen, Mi{H), divided by the white dwarf's cooling age, tcooi,i- 



MAH) 



tcooLi 



(3) 



As explained in I Jural (|201ll ). we expect the loss of accreted hydrogen by the action of a stellar wind 
to be negligible. To determine Mi{H), we have used reported values of [H]/[He] and models for 
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the mas s of the s t ar's c onvective zone as a function of stellar gravity and effective temperature. 
Because iKoester I ((20091) only presented results for one value of white dwarf gravity and hydrogen 
composition, D. Koester kindly provided a grid of models for DB stars with a sufficient range of 
temperature, mass and [H]/[He] to enable us to interpolate for the mass of the convective zone for 
each star individually in Table 1. The estimates of Mi{H) are strongly sensitive to the star's total 
mass. For example, for DBs with an effective temperature of 13,000 K, the mass of the convective 
zone is nearly a factor of 200 smaller in a star of 1.1 Mq compared to a star of 0.6 M ^ ti (IDufour et al. 



2010l ). In most cases, we take [H]/[He] from the optical study of iBergeron et al. I (|201lk but for a 
few stars hotter than 20,000 K, we use upper bounds from published ultraviolet observations which 
are appreciably more sensitive to the hydrogen abundance. All our e stimates for MjjH) ar e upper 



bounds because the white dwarf could possess primordial hydrogen (jBergeron et al. 



2011 



As shown in F igure 1, for the seven stars in Table 1 for which both IBergeron et al. I ()201ll ) and 
Voss et al. I (j2007l ) derived masses spectroscopically, instead of just adopting an assumed mean. 



the values of Mi{H) are systematically lower from our use of the IBergeron et al. I (|201ll ) study 
compared to those in I Voss et al. I (|2007l ) by an average factor of 4 . This substantial disc repancy is 
the result of four factors all operating in the same direction. First, IBergeron et al. I (|201ll ) typically 
found larger stellar masses. The typical 15% increase in derived mass means that the mass of the 
convective zone is smaller by approximately a factor of 2. Second, stars with a larger mass take a 
longer time to cool to the observed effective temperature; this effect might contribute a factor of 
1.2 to the estimate of MjjH). T hird, because the star's gravity is estimated to be higher in the 
analysis of lBergeron et al. I (j201ll ). the ratio of the hydrog en to helium may be found to be as much 



as a factor of 2 lower than derived by I Voss et al. I (|2007l ) . Finally, there is a different theoretical 



treatment of convection such that the valu e s of t he mass of the convective zone are typ i cally a 



factor of 1.5 lower in the models of iKoester I (j2009) compared to those us ed bv lVoss et al. I (|2007l ) 



The spectroscopically-derived masses derived by IBergeron et al. I (j201ll ) usually agree very well 
with those derived from trigonometric parallaxes, and therefore it seems likely that the values of 
Mi{H) derived from their analysis are realistic. It must be recognized, however, that there are 
uncertainties. 



A few stars in Table 1 have remarkably low upper limit s to MjjH ). These stars can be used 
to place a bound on the space density of interstellar comets ([Jural l201ll ). 



2.3. Heavy Atom Accretion Rates 

We wish to compare the hydrogen and heavy atom accretion rates. For the i'th star and the 
j'th element, the time- averaged mass accretion rate from the disk onto the star, Mi{Zj), we take 



(4) 
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where Mi(Zj) is the inferred mass of element Zj in the star's mixing zone where the gravitational 
settling time is ti{Zj). Using Ztot to denote the mass of all the heavy elements, we define for each 
star the time-averaged total heavy element accretion rate, Mi{Ztot), as: 



MiiZ, 



tot 



(5) 



In many cases, only calcium is measured because this element is the most easily to detect in optical 



spectr a; the total heavy element accretion rate is extrapolated as described bv lZuckerman et al. 



(|201Cll ). The values of Mj i Ztnt) are relatively insensitive to the white dwarf mass because, as 



illustrated in iKlein et al. I (120101 ). the settling time scales approximately with the mass of the 



convective zone. 



Three sources of hydrog en in DBs have been proposed: inters tellar (IVoss et al. 1120071 ). cir- 



cumstellar (jjura fc: Xu Il2010l ) and primordial (jBerg eron et al. II2OI1I ). Because the heavy elements 
settle while the hydrogen does not, the absence of a correlation in Figure 2 is not a decisive test of 
any particular model for the source of the DBs' hydrogen. 



2.4. Summed Accretion Rates 

Although there are nearby DBs that are missed, we assume that in identifying the stars, there 
is no bias related to their external-pollution. If so, then we can treat the set of stars as a well 
defined ensemble and therefore compare the average hydrogen and heavy atom accretion with each 
other in a meaningful manner. We define Msum{H) as: 

Msum{H) = Y^miH) (6) 

i 

For comparison, the summed heavy atom accretion rate, MsumiZtot), is: 

MsumiZtot) = Yl MiiZtot) (7) 



With estimates in Table 1 of Mi[Ztot) for 2/3 of the stars, and ignoring any contribution from 
the remaining 1/3 of the stars, Msum{Ztot) is 1.6 x 10"'^'^ g s^-*^. Nearly all stars in Table 1 have 
been examined for hydrogen, and Msum{H) is <1.4 x 10^ g s^^, approximately a factor of 1000 
smaller than the rate for heavy atoms. Because some currently-unknown fraction of this hydrogen 
is derived from the interstellar medium or is primordial, we conclude that in the entire ensemble 
less than 1% of the accreted mass is carried in water. This bound is only strengthened if there 
is appreciable accretion of heavy atoms in those stars which have not been examined with high 
spectral sensitivity. Because we examined 57 stars, the upper bound to the average accretion rate 
for an individual star is 2.5 x 10^ g s~^. 
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For both heavy elements and hydrogen, the asummed accretion is dominated by a few stars. 
If we remove the top three accretors in each category, then Msum{Ztot) and Msum[H) are <2.6 x 
10^ g s~^ and 4.4 x 10^ g s~^, respectively. In this analysis, the ratio of the two accretion rates 
is approximately a factor of 600. Considered as an aggregate and interpreting the data with the 
best available parameters, we find that extrasolar asteroids accreted onto the warm DB stars are 
appreciably drier than CI chondrites. 

For the four stars in Table 1 for which comprehensive abundance studies have been performed, 
we show in Figure 3 a plot of the fraction of the total mass of each of the major individual 
constituents - O, Mg, Si, Ca and Fe - and the upper bound to the fraction of the toal mass that 
could have been water. In this Figure, we assume the steady state approximation given by Equation 
(4). It can be seen that there is a relatively small scatter in most of the mass fractions - especially 
oxygen, consistent with our treatment of the extrasolar asteroids orbiting separate stars as one 
population. We also see that the upper limit derived from the analysis of the aggregate population 
is significantly lower than the limit for any individual star. The individual bounds plotted in Figure 
3 are derived by using the lower of the two values for the water fraction in the parent body. The first 
value which is usually the stronger contraint assumes that all the oxygen in excess of that which 
could be bound into MgO, Si02 and Fe203 was contained in water. The second value assumes all 
the atmospheric hydrogen was bound into water. As explained in Section 3, the evidence that CD 
61 has accreted ice-rich material is ambiguous, and we therefore plot the water fraction as an upper 
limit rather than measured quantity. 

We see in Figure 3 that the extrasolar asteroids a re in aggregate drier even th an CV chon- 



drites - the most refractory-rich of common meteorites (|Wasson &: Kallemevn Ill988l ). We also see 
that commonly the extrasolar asteroids also are more rich in the refractory element Ca than the 
CV chondrites. We can understand the dryness of extrasolar asteroids as being a result of their 
having formed interior to the snow line. We do not yet have a good model to explain the relative 
enhancement of Ca that is measured in these systems. 



3. MODELING INDIVIDUAL ACCRETION EVENTS 

Having established statistically that ice is low in abundance within our sample of extrasolar 
asteroids, we now revisit model interpretations of elemental abundances in individual stars. Because 
different elements settle at different rates, there is not necessarily a simple proportionality between 
the rel ative abundance s in the atmosphere of the white dwarf and the fractional masses in the parent 



body. iKoester I (|2009l ) described three regimes for an accretion event. Initially, after an asteroid 
is tidally-disrupted into a disk, there is a build-up of heavy elements in the outer convective zone 
of the star. In this regime, the relative abundances in the stellar atmosphere directly scale as the 
abundances in the parent body. However, once the elapsed time since the onset of accretion becomes 
comparable to the gravitational settling time, the system enters a second phase. In this situation, 
it has been assumed that the gravitational settling rate balances the accretion rate and an exact 
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steady state is established. Finally, the disk material is completely accreted and the system enters 
a third phase when the atmospheric pollution decays with the elements with longer settling times 
lingering in the atmosphere and nominally appearing overabundant. During the first and second 
phases, the white dwarf possesses a cricumstellar disk. In the third phase, the disk is dissipated. 
If the disk i s dusty, then we can detect an infrared excess. If, however, the disk is largely gaseous 
(jJura II2OO8I ) then it may not be observationally evident. 



From the modeling point-of-view, the second phase is the most uncertain because material is 
both being added and lost from the mixing zone. In the first phase, material is only being added and 
in the third phase, material is only being lost. We therefore reconsider models for the second phase 
of an accretion event and reconsider the assumption of an exact steady state. T he physics of the 
accretion from the disk onto the star is not fully understood; it might be variable. iRafikov I ()2011al ) 
has argued that Poynting-Robertson drag on the disk is important, but clearly this is only part o f 
the story and the accretion rate may vary significantly (jRafikov Il2011bl : iBelvaev fc: Rafikov II2OIII ) 
during a disk's evolution . Observations of the time-variat i on and profiles of emission lines from 



the circumstellar matter (IGaensicke et al. 



2006 



20071 . l2008l : lMelis et al. II2OIOI ) raise the possibility 



that there are significant changes in the accretion rate. We now consider a quasi steady-state where 
there is a disk, but, nevertheless, the accretion rate varies in time and the system is not in an exact 
steady state. 



Dupuis et al. I () 19931 ) computed models for accretion of heavy atoms from the interstellar 



medium onto DB white dwarfs. Because they invoked both high rates of accretion from clouds and 
low rates of accretion from the intercloud medium, their calculations explored the observational 
consequences of time-varying accretion and showed that the photospheric ratio of Mg to Ca could 
be dramatically enhanced because Mg gravitationally settles more slowly. While recent evidence 
demonstrates that the heavy atoms in DBs mainly originate from asteroidal parent bodies, their 
finding that the star's atmospheric abundances are sensi tive to the time - histor y of the entire system 
remains valid. Our models extend the previous work of lDupuis et al. I (j 19931 ). 



3.1. Quasi-steady Accretion Model 

In the outer mixing z one of a white dwarf, the balance between accretion and settling is 



governed by the expression (IKoester II2009I ): 

M(7A 

+ Mpb{Zj) (8) 



dt tj 



where Mps^Zj) denotes the accretion rate from the circumstellar disk and is a measure of the 
composition of the tidally-disrupted parent bod}0. Performing an integration on both sides and 



^For notational convenience, in this Section, we denote the mass in the mixing zone in each star as M,{Zj) instead 
of Mi{Zj) as in Section 2. 
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assuming that M^,[Zj) = at t = 0, the solution to Equation (8) is: 





t/h I e^'/i^ MpB{Zj)dt' 



(9) 



Inverting the integral in Equation (9) is required to determine the parent body composition. One 
approach is to Fourier analyze the accretion rate. Because of the large unknowns, for simplicity, 
we consider just the illustrative special case of a single Fourier term: 



MpB{Zj) = MpB{Zj) (1 + r sina;t') 



(10) 



Here r and oj are free parameters that characterize the amplitude and frequency of the fluctuations 
in the accretion rate of the j'th element whose time-average value is MpB^Zj). In addition to 
accretion rate variations intrinsic to the disk, there may b e multiple d isruptions of parent bodies 
that also may lead to time variations of the accretion rate (jJura II2008I ). In all cases, we require r 
< 1. 

For convenience, define the dimensionless term: 



Bj = tjuj 



Then: 



M,{Zj) = MpB{Zj) tj 



+ 



sinwt — Bj cos cot + Bje 



(11) 



It is instructive to consider some limiting case solutions to Equation (12). When t « tj, then: 



M,{Zj) ^ MpB{Zj 



(13) 



In this circumstance, the element's mass in the stell ar atmosphere d irectly corresponds to the 
accretion rate; the initial build-up phase discussed bv iKoester I (|2009l ). When t » tj, then the 
solution to Equation (12) is: 



M,{Zj) « MpBiZj) tj 1 + 



1 + B] 



[sin ut — Bj cos uj t] 



(14) 



In the exact steady state where the accretion rate is constant and r = 0, then the abundance of 
an element in the parent body are just given by the abundances in the photosphere divided by the 
settling time and from Equation (14) we recover Equation (4). However, if r > 0, the relationship 
between the photospheric concentration of an element and its relative abundance in the parent 
body is not the simple proportionality of Equation (4). 

By examination of Equation (14), we can see that the quasi steady-state approximation usually 
agrees with the exact steady state model by better than a factor of 2. That is if -Bj >> 1, then clearly 
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M*(Zj) is essentially constant and the limit of Equation (4) is reached. In this case, variations in 
the accretion time are short compared to the settling time and the pollution mass in the mixing 
zone is constant. If Bj « 1, then M^,{Zj) might vary substantially, but the variations in the 
accretion rate are so slow that the system's behavior is very similar to the exact steady state. The 
only situation when there is a significant difference between the quasi-steady state and the exact 
steady state is when Bj 1. In this case, the settling and accretion times are comparable and the 
interplay between the two factors can be complex, each element behaving differently according to 
its specific value of Bj . 



3.2. Was the Parent Body Accreting onto GD 61 Ice-Rich? 

We now use the model of quasi-steady accretion to a ddress the question o f whether the parent 



body accreted onto GD 61 was ice-rich as proposed bv iFarihi et al. I (j2011al ). We see in Figure 
3 that both G241-6 and GD 61 have relatively high fractions of oxygen in their contaminants, 
and therefore both systems are candidates for the accreted parent body having had a substantial 
amount of water. However, G241-6 has little photospheric hydr ogen, and therefore less than 10% 



of the mass fraction of t he pollution is wa ter (iKlein et al. Il201ll ). Because G241-6 does not have a 



circumstellar dust disk (IXu k. Jura Il201ll ). plausibly, it is in a late phase of an accretion event and 



the oxygen is especially abundant because it lingers longer in the outer settling zone. In contrast. 



(Farihi et al. 


2011a^ 


Farihi et al. i 


2011al 



assume that the GD 61 system is in an exact steady state. If so, then by mass 
there is about 30% more oxygen than can be locked into oxide minerals and this implies that the 
parent body contained a substantial amount of ice. This oxygen excess is sufficiently small that the 
analysis is sensitive to whether the system is in an exact steady state or only a quasi steady state. 

One difficulty with using the exact steady state model to interpret the data for GD 61 is that 
iron displays a relatively low abundance. In an exact steady state model, as shown in Figure 3, it 
is inferred to have an usually low fractional abundance compared to other polluted white dwarfs. 
One possible way to explain simultaneously both the low iron and the high oxyge n fractional 



Klein et al. 



2011 



aburidances is that the system simply is in a late phase where only settling occurs (jJura fc: Xu 
2010 



. However, this hypothesis is inc onsistent with the pr esence of a dust disk 
and therefore the likelihood that accretion is ongoing. iFarihi et al. I (|2011al ) suggested that the 
accretion onto GD 61 is iron-poor because the circumstellar disk is composed of the outer portion 
of a planetesimal where most of the iron was concentrated into a central core. 



GD 61 is deficient in carbon relativ e to the Sun by ab out a factor of 1000 (i Desharnais et al 



20081 : IFarihi et al. Il2011al ). According to lLee et al. I (j2010l ). in planet-forming disks, carbon should 
be treated as a volatile which therefore explains w hy this element is commonly observed to have a 
low abundance in extrasolar asteroids (I Jura II2006I ). If, in fact, the parent body accreted onto GD 
61 was ice-rich, it would be a puzzle to understand why it was simultaneously carbon-poor. 
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Here, we suggest that a quasi-steady state model may explain the data without any requirement 
that the accreted parent body contained ice. Therefore, our input parameters are chosen to find the 
largest effect of the quasi static model; they are not necessarily physically realistic. Consequently, 
using Equation (14), we compute a model with uj = 0.8 x 10~^ yr~^ a, nd r = 0.95. Because the 
gravitational settling times typically are ~10^ yr ( Farihi et al. Il2011al ). this choice of co means 
that the values of Bj are nearly 1. The relative mass fractions of O, Mg, Si and Fe are taken to 
equal 0.40, 0.18, 0.2 and 0.22; these relative abundances are close to the values for bulk Earth 



([Allegre et al. Il200ll ) except for iron which is about 50 % too low. However, because the error in 



the photospheric abundance of Fe in GD 61 is 0.2 dex (jFarihi et al. Il2011al ). the true fractional 



abundance of this element in the parent body is not required to be anomalously low compared to 
bulk Earth. 

We display in Figure 4 the results of our model calculation for the fractional abundances in 
the photosphere at different times. We see that there is a phase when the model matches the data. 
It therefore seems that the argument that there must be water in the parent body accreted onto 
GD 61 is provisional. 



DISCUSSION 



Hydrogen and oxygen are sufficiently cosmically abundant that in a disk where gas condenses 
into solids at least half of the total mass could be ice. If, however, the temperature is high enough, 
only relatively refractory materials enter the solid phase and eventually are incorporated into plan- 
etesimals. Because the temperature in a disk decreases radially from the central star, it is likely 
that ice-rich objects form in cold, distant regions. We find that water is less than 1% of the mass 
of the ensemble of extrasolar asteroids accreted onto DB white dwarfs, and therefore these plan- 
etesimals likely formed interior to the snow line, the theoretical boundary between the ice-forming 
and the ice- free zones. 

We have argued that oxygen is not sufficiently enhanced in GD 61's heavy element contami- 
nation to demonstrate convincingly that there was water within the parent body. In the future, a 
polluted white dwarf might be identified where the oxygen is sufficiently more abundant than the 
other species that the argument for water could be more convincing. In GD 40's photosphere, the 
numbe r of O atoms is a f actor of 2.3 times greater than the sum of the number of Mg, Si and Fe 
atoms (iKlein et al. Il2010l ). In this star, there is no evidence for any ice in the accreted parent body. 
In GD 61's photosphere, the num ber of O atonis is a f actor of 2.9 times greater than the sum of the 
number of Si, Mg and Fe atoms (jFarihi et al. Il2011al ). In this star, there may have been water in 
the parent body, but the case is uncert ain. In the Suri , the number of O atoms is 4.9 times greater 
than the sum of Mg, Si and Fe atoms (jLodders II2003I ). A planetesimal that formed in a very cold 
environment could have a correspondingly large fraction of oxygen. If this object is accreted onto 
a white dwarf, its very high oxygen fraction could be measured and the case would be good for an 
ice-rich parent body. However, if the white dwarf is a DB, then the settling times are sufficiently 
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long that unless there was also evidence for ongoing accretion, it might be difficult to exclude the 
possibility that the enhancement of oxygen is simply the consequence of its lingering longer in the 
mixing zone. 

We could hope to use warm DA white dwarfs to study relative oxygen abundances because 
the settling tim es are shorter than a year and the system almost certainly is in a steady state 
(jKoester II2009I ). Unfortunately, the current situation is murky. GALEX J 1931 is an exceptionally 
heavily polluted DA where all the major element contaminates are measured, and the number of 
O atoms i n the photosphere is somewhere between 1.0 and 1.8 times the sum of the Mg, Si and 
Fe atoms (|Vennes et al. Iboiol . [2011al : iMelis et al. Iboilh . Th e two groups d o not agree very well 
in their determinat ion of the iro i i abun dance, perhaps because I Vennes et al. I (|20ld ) only used one 
weak Fe line while iMelis et al. I (j201ll ) used several iron lines. More importantly for our purposes 
here, correcting for settling and consequent stratification through the atmosphere, the num ber of 



O atoms conipared to the sum of Mg, Si and Fe atoms is either ~11 (jVennes et al. II2OIOI ) or ~1 



( Melis et al. 



2011 



In the first case, likely the parent body was ice-rich; in the second case likely 
it was ice-poor. Because we have found that ice-rich parent bodies are not ubiquitous; observations 
of additional DA stars may help resolve this disagreement. The analysis of DB stars does not suffer 
this particular ambiguity. 

In our toy model for external pollution, we have assumed that the ultimately-accreted hydrogen 
is chemically bound to oxygen in the form of water. Another possibility is that this hydrogen could 
be chemically bound to carbo n in the form of hydrocarbons. However, carbon is highly volatile in 
planet-forming environments ( Lee et al. 20ld) and deficient in tho se relatively few polluted white 



dwarfs where it has been studied (Jura 



20061 : iFarihi et al. II2009I ). Therefore, H2O is the most 



likely carrier of hydrogen in extrasolar minor pl anets. This water could be contaiii ed in hydrated 
silicates such as serpentine rather than pure ice (jRivkin et al. 1120021 : IWasson II2008I ). 



A large source of uncertainty in our estimates to the upper bounds to the hydrogen accretion 
rates is the accuracy of the adopted white dwarf masses. This difficulty should be eased with 
the Gaia mission which will determine trigono metric parallaxes of a ll nearby white dwarfs, greatly 



extending the work of the Hipparcos satellite (jVauclair et al. 1119971 ). 



5. CONCLUSIONS 

We argue that the aggregated external pollution from parent bodies accreted onto DB white 
dwarfs is less than 1% by mass composed of water. Consequently, on average, the asteroids that 
orbited main-sequence stars of typically 2.5 M0 are drier than CI chondrites by at least a factor 
of 20. In contrast, we show that for an individual star such as CD 61, having a large fraction 
of oxygen in the polluted material does not necessarily imply accretion from a water-rich parent 
body. Despite the ambiguity associated with studying any single white dwarf, our analysis of the 
aggregate population of polluted DBs provides indirect observational support that snow lines in 
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planet-forming disks are common. 

We thank D. Koester for providing a grid of models for the masses of white dwarfs and for 
useful email exchanges. We thank the anonymous referee for thoughtful and constructive comments. 
This work has been partly supported by the National Science Foundation. 
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Table 1 - DBs Within 80 pc of the Sun 
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Table 1 - Continued 
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The entries in the column for M{H) may be blank (no measurement reported in the literature), < 
(no hydrogen detected, but an upper limit is reported) or < (hydrogen detected, but the rate of ac- 
cretion is an upper limit because there could be some primordial hydrogen). We enter "<" in the log 
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M{Ztot) column for those stars examined in the SP Y survey that do n ot display any calcium absorp- 
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Fig. 1. — Comparison of time-averaged upp er bounds to the h ydrog en accretion rate s , den oted 



M{H) in the text, using the parameters in IVoss et al. I (120071 ) and iBergeron et al. I (120111 ) for 



the seven stars in Table 1 for which both studies derived masses spectroscopically. The dotted line 
displays the locus of points where the two rates would agree. As discus sed in the text, the difference 
between the two studies is largely the result of IBergeron et al. I (|201ll ) deriving higher white dwarf 
masses. The errors are only those associated with the atmospheric abundance d etermination s ; they 
do not reflect uncertainties in the stellar mass or effective temperature. While IVoss et al. I (|2007l ) 



found that the uncertainties in their derived hydrogen abundances might be as low as 6%, our 
plotted error bars reflect our assumption that realistically, the errors may be as large as 0.1 dex. 
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Fig. 2. — Comparison of the upper bound to the time-averag ed hydrog en accretion rates, denoted 
M{H) in the text, with heavy atom accretion rates, denoted M{Ztot) in the text, for the stars in 
Table 1 where both quantities are determined. For the stars where only the calcium abundance 

is reported, we assume an overall uncertainty of a factor of two in the total heavy-atom accretion 
rate. For the stars with detailed abundance analysis, the plotted error bars are taken from the 
papers where the data are reported. No correlation between hydrogen accretion and heavy atom 
accretion is evident. 
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Fig. 3. — The relative fraction of the total mass of the accreted parent body carried by individual 
constituents for the pohuted white dwarfs in Table 1 where all dominant elements - O, Mg, Si, 
Ca and Fe - have been measured. The upper limit labeled with "Sum" is the bound placed 
on the fractional aggregate water content of extrasol ar asteroids in Secti on 2. The abundances 



assigned to the parent bodies accreted onto GD 61 (IFarihi et al. 112011a), GD 40 (jKlein et al 



2OIOI), G241-6 (IZuckerman et al. Il20ld ) and HS 2253+8023 (|Klein et al. 



2011) are all corrected 



from the atmospheric abundances by assuming a steady state and therefore using Equation (4). 
The water fractions for the white dwarfs are all presented as upper limits; the case of GD 61 is 
discussed in detail in Section 3.2. The water fractions for solar system CI and CV chondrites are 
shown for comparison; they are computed by assuming all the hydrogen in the meteorites is carried 
in water in the form of hydrated minerals. 
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Fig. 4. — Predicted relative mass fractions of individual heavy elements from the accreted parent 
body within the photosphere of GD 61 in the quasi-steady state model given by Equation (13) with 
the parameters prov ided in the text. The dotted horizontal lines display the values measured by 



Farihi et al. I (|2011al ) and the solid curves the prediction as a function of time scaled to dimensionless 



units. Even though there is no ice in the model parent body, at ojt ~ 4.9, the data are well matched 
by the model. The computed fraction of the contamination that is oxygen is always greater than 
the assumed value in the parent body of 0.4 because oxygen settles relatively slowly compared to 
the heavier elements. 



